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Astrophysical models for the high-energy emission of blazars are reviewed. Blazars ejecting 
relativistic radio jets at small angles to the line-of-sight are the only type of active galactic 
nuclei (AGN) discovered above 100 MeV. The 7-rays apparently originate in the jets which 
explains the absence of 7-ray pair creation attenuation. The bulk Lorentz factor of the radio 
jets is much smaller than the Lorentz factors of the emitting particles requiring in situ particle 

■ acceleration such as Fermi acceleration at shocks. The mounting evidence for a correlation 
■^j- ' between the optical and 7-ray emission argues for the same accelerated electrons emitting 

the polarized optical synchrotron photons to be responsible for the high-energy emission. 
Alternatively, the 7-rays could be due to electrons of a secondary origin related to the energy 
losses of protons accelerated at the same shocks as the synchrotron emitting electrons. In this 
case blazars would produce an observable flux of high-energy neutrinos. Unified schemes for 
AGN predict a circum-nuclear warm dust torus attenuating 7-rays above ~ 300 GeV emitted 
from within a central sphere of radius ~ 2 x 10 4 rs which rules out external Compton scattering 

■ models as the origin of the TeV 7-rays from Mrk421 and Mrk501. 

a: 
o 

5h ' 1 Introduction 

The EGRET discovery of 7-ray sources at high galactic latitudes and their identification with 
blazars (~ flat-spectrum radio-loud AGN) has prompted a large number of authors to think 
about possible explanations for the phenomenon. As of 1996, I have counted almost one model 
or one model ramification per source scanning through NASA's Astrophysics Data System. The 
considerable theoretical interest documents the importance of the discovery and at the same 
time demonstrates the on-going struggle for a viable explanation. Of course, not all of the 
papers can be discussed in this review, and I have decided to critically discuss only the basic 
model assumptions and their consequences. The assumptions that (i) blazars represent accreting 
supermassive black holes ejecting bipolar jets along the angular momentum axis of the inner 
accretion disk and that (ii) the 7-rays originate in the jets are common to all the models. This 
consensus is based mainly on the facts that EGRET has observed 7-rays only from AGN with 
jets and that jets provide a natural resolution of the so-called compactness problem: For 7-rays 
of energy m e c 2 the pair creation optical depth 7 + 7 — > e + + e~ of an Eddington accreting black 
hole is the universal constant 

„ L cdd 3a T = 3 m p _ ^ 
77 47rr s m e c 3 16 32 m e w 

Eddington accretion onto a supermassive black hole is the only mechanism efficient enough 
to explain the high luminosities and short variability time scales of blazars assuming isotropic 
emission. However, th^ large optical depth Eq.(l) is inconsistent with the observed optically 
thin power law spectrs-0 and therefore the emission must be anisotropic. 
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Figure 1: Coarse overview of models. See Sect. 1.1 for details. 

The accreting black hole picture for blazars is based on unified schemes for AGN which are 
discussed extensively in the contribution by Padovanca. As a general cautious remark I want 
to emphasize that the predictions of the unified schemes for the density of matter and radiation 
in the medium surrounding the jets are subject to change in the course of further astronomical 
studies and should therefore be regarded as preliminary. Furthermore, the unified schemes leave 
it essentially open whether the jets form as a hydromagnetic e~p wind from the relativistic 
inner accretion diskUJ or as a Poynting flux driven wind extracting rotational energy from the 
central black holJ]. The merging history of ellipticals as th e t ypical host galaxies of radio-loud 
AGN lends some support for the rotating black hole pictured, although jets are also associated 
with young stars which do not contain black holes. It is hoped that the 7-ray properties of 
blazars probing the physical conditions in the jets ameliorate theory in this important issue. 
However, baryon pollution^ could also wash out the difference between pair winds and e~p disk 
winds. 

1.1 Relativistic particles, cold and hot beams, internal and external target photons 

By analogy with the local cosmic rays, the relativistic particles responsible for the 7-rays are 
expected to be electrons and protons. In the tenuous plasma of the jets, electrons (positrons) 
produce 7-rays mainly by inverse-Compton scattering or synchrotron emission. For protons 
these processes are generally unimportant owing to their larger mass. The pion production 
on matter, which is important as a 7-ray production process in the Galaxy, can be neglected 
due to the low plasma density in radio jets. However, the photo-production of secondary pions 
and pairs in collisions with ambient low-energy photons can be important. In spite of the 
hadronic cross-section <j p7 = 5 10 -28 cm 2 being much smaller than the Thomson cross-section 
o"t = 6.7 10~ 25 cm 2 , the proton cooling time scale can be as short as the electron Thomson 
cooling time scale owing to the much larger inelasticity of hadronic processes. The proton and 
electron energy losses (i = e,p) are given by 



where the factors / e = 1 + u 7 /mb and / p = 1 + 500u 7 /mb take into account the relative contri- 
bution of synchrotron, inverse Compton, Bethe-Heitler pair production, and pion production. 
Balancing with energy gains due to Fermi acceleration 

§=K«?B (3, 
where K < 0.4 depends on the diffusion coefficients one obtains the maximum Lorentz factors 

7p,max = 710 10 iO J Er5/-5 ( 4 ) 

for protons and 

7e,max=4l0 7 J ?d£r5/ p ^ ( 5 ) 

for electrons. Adopting K = 0.1, B = 1 G, and u 7 = ur this corresponds to maximum energies 
-^p,max = 10 6 TeV and E ejJa3X = 4 TeV. The assumption of K = 0.1 is certainly an overestimate 
for electrons which, owing to their much smaller Larmor radius, are in resonance with much 
smaller scales of the expected inverse-power-law plasma turbulence spectrum than protonfi. 
The ratio of cooling time scales for electrons and protons is given by 
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(6) 



from which one can readily approximate the proton induced luminosity in terms of the electronic 
inverse-Compton and synchrotron luminosity 

L p ^L e = ^(^y^L e (7) 

M e t p M e \m p J 7e / e 

where u p and u e denote the energy densities in protons and electrons, respectively. To estimate 
the ratio of Lorentz factors relevant for the luminosities it is important to realize that the 
synchrotron luminosity is mainly due to electrons with t e < tdyn where idyn = r sh/c denotes 
the dynamical time scale associated with the expansion of the shock front. The luminosity is 
not dominated by emission from the most energetic electrons, but increases only logarithmically 
above t e = idyn corresponding to Lorentz factors 7 e ~ 10 2 (infrared break). With 7 Pimax ~ 10 9 
one obtains from Eq.(7) a maximum proton- induced luminosity of the order of 

ip.max ~ —L c (8) 
M e 

Photoproduction by ultra-relativistic protons in a relativistic jet has been considered by a num- 
ber of authors assuming different locations of the 7-ray emitting zone, either very close to 
the central black hokBH, at a distance of the order of the broad emission line region radiu@, 
or beyond the broad line region where the jet is dominated by its own synchrotron radiation 
Interestingly, the proton maximum energies and luminosities required to explain the 
blazar 7-rays conspire to produce the observed flux of ultra-high energy cosmic rayS 

The bulk flow of relativistic particles can be distinguished as being either a cold beam or a 
hot beam. In the former case, the particles have zero kinetic energies in the frame co-moving 
with the bulk flow. The TeV 7-ray observations would then require cold beams with particle 
energies of at least several Te\S. In the latter case, the bulk flow need only be moderately 
relativistic to accommodate for the superluminal motion and the absence of pair attenuation, 
whereas the co-moving frame energies are ultra-relativistic. Superluminal motion of radio knots 
suggests bulk Lorentz factors r ~ 10 (at the ~ 0.1 pc scale) and observations of the synchrotron 



turnovers indicate co-moving frame Lorentz factors of up to several tens of_GeV. Cold beams 
could be generated through ordered electric fields c. 
fields arising at standing shock waves! Cold beamed 



the central objedEffl or electric drift 
are expected to slow down reaching 



bulk Lorentz factors of T ~ 10 owing to Compton dragHEI. Moreover, the same target radiation 
field which is assumed to be up-scattered to 7-ray energies would attenuate the 7-rays by pair 
production on a comparable time scale above photon energies of m e c 2 . An interesting ideaa is 
that cold sub-beams embedded in a larger scale jet could be produced at shocks in the winds of 
stars orbiting through the jet plasma. To explain the variability in Mrk421, one would, however, 
need an unexpected high density of mass-loosing stars in the central stellar cluster. Since the 
Compton-drag on cold beams is difficult to avoid and the bulk Lorentz factors inferred from radio 
observations are much smaller than the Lorentz factors of the emitting particles, one needs in 
situ particle acceleration further downstream the jet as assumed in the hot beam models. Hot 
beam models with only one zone of emission in the jet (a 'blob') and external photons as the 
dominant target for inverse-Compton scattering (external Compton models: EC) have been 
published assuming that the blob is close to the nucleullil at a distance of (10 2 — 10 3 )rs or at 



the distance of the broad emission line region e3 at (10 3 — 10 4 )rs- In the hadronic variant of 
this idea?!, protons are assumed as the primaries initiating syn chrot ron cascades by scattering 
external photons (E-PIC). In synchrotron-self-Compton modeli?iHi3'L3, the 7-rays originate in 
the jet where it is dominated by the internal synchrotron photons, i.e. presumably at distances 
of > lQlTs- Again, there is a hadronic analogue, the proton-initiated synchrotron cascade 
modePIr 1 ! (PIC, denoted as S-PIC in Fig. 4). The closer the 7-ray emitting zone to the center, 
the lower the 7- ray energy which can be emitted without attenuation. The inhomogeneous 
EC modelsO'EKa consider the emission of a Poynting flux/pair jet from close to the nucleus to 
the edge of the broad emission line region, the effect of this energy-dependent 7-photosphere is 
that with increasing photon energy, the 7-rays come from increasingly distant parts of the jet. 
Recently, it has been pointed oucil that the high density of pairs required in these models at 
r < 10 2 rs would lead to a bump in X-rays from inverse-Compton scattering of disk photons. 
Instead of a bump, blazars show a deep gap in X-rays. As can be seen from Fig. 1, all possible 
combinations of the basic assumptions have been investigated. In Sect. 2, the model predictions 
are confronted with observations. Section 3 concludes with a critical discussion of the models 
for high-energy emission from blazars. 



2 Key observational results and model predictions 

2.1 Fraction of 7 -ray sources among blazars 

A fraction of ~ 10% of all blazars have been detecteoH in the EGRET pass-band 100 MeV - 
10 GeV. The blazar fraction refers to those flat-spectrum radio sources which would have been 
above the EGRET flux sensitivity if the ratio between their 7-ray and 5 Ghz radio fluxes were 
at least as large as the lowest value for the EGRET detected blazars. An obvious explanation 
of the 10% would be flux variability which is, however, inconsistent with the 30% duty cycle 
of the strong EGRET detections. Intermittent behavior is suggested by interferometric radio 
observations which show months or years between the ejection of new radio knots, possibly 
associated with periods of 7-ray flaring and quiescence. Nevertheless, some quasi-stationary 
emission models make definite predictions about the 7-ray blazars representing a true subset of 
all blazars. EC models predict a stronger Doppler boosting of the 7-ray flux oc <5 4 + 2a compared 
with the radio-to-optical synchrotron flux which is boosted only by the factor 5 3+a . With a = 1, 
the prominence of the 7-ray component increases oc 5 2 implying that the most pronounced 
7-ray blazars must be the ones with the largest Doppler factors (smallest viewing angles or 
largest bulk Lorentz factors). This prediction can be tested experimentally with interferometric 



radio observations. There are claims that, indeed, the strong EGRET sources have atypically 
high values of the Doppler factor Jl, whereas there the Doppler factors of the EGRET sources 
among the S5 radio sources are normal for blazarsJ. In the PIC model, the ratio between 
7-ray and radio-to-optical luminosity depends on the proton-to-electron ratio and on the ratio 
of cooling time scales, see Eq. (5). There is no reason to believe that either one should be 
constant, it is rather plausible that the numbers vary with shock obliquity, Lorentz factor, and 
jet magnetization. Proton acceleration is likely to be limited due to the particle gyroradius tl 
exceeding the shock radius of curvature r s h rather than due to energy losses. If this is indeed 
the case, one obtains tl oc r s h oc rj at the maximum energy with rj denoting the jet radius. 
From eq.(7), the ratio of proton induced 7-ray and electronic synchro/Compton luminosity is 
then given by L p /L c oc 7 Pi max oc r^B. Jet formation models give the scaling rj oc Mbh and 

1 /2 —1/2 

B ^cdd-^bh (from ub = w c dd)- The parameter Z e dd = -^/Asdd denotes the luminosity in 

1/2 1 /2 

units of the Eddington luminosity. Hence one obtains L p /L e oc M b £ Z odd . Thus, in the PIC 
model quasars, which are intrinsically more luminous than BL Lacs and therefore require larger 
Mbh's or / e dd' s ) are expected to be relatively more luminous in 7-rays than BL Lacs which is 
consistent with observations. The non-detected blazars would correspond to jets associated with 
sub-Eddington accretion flows. 

2.2 Amplitude of 7 -ray /radio luminosity ratio 

The ratios between the observed 7-ray and radio luminosities range roughly between 10 2 and 10 4 
(Fig. 2). If 7-ray blazars have a true distribution of 7/radio amplitudes larger than 100, so that 
only the 10% with the largest amplitudes show up above the EGRET sensitivity threshold, the 
observed blazar fraction could easily be explained. In the EC model, the amplitude distribution 
is triggered by the jet Doppler factor and with 5 = 1 — 10 the model predicts amplitudes between 
5 2 = 1 — 100. Maximum Lorentz factors somewhat larger than 10 would therefore be sufficient 
to explain the width of the amplitude distribution. As shown above, the PIC model requires 
that the product of black hole mass and Eddington ratio spans a range of more than 10 4 which is 
also necessary to explain the range of blazar bolometric luminosities. Whereas the EC and PIC 
models entail satisfactory explanations of the blazar fraction and the amplitude distribution, 
the SSC models are inconsistent with the observations of a few high-luminosity quasars such as 
3C279 and PKS 0528-134 in which the 7-ray luminosity exceeds the radio-to-optical synchrotron 
luminosity by factors 10 — 100. This follows from re-writing the ratio of the singly scattered 
Compton luminosity to the synchrotron luminosity in terms of the Thomson optical deptlc3 
which yields Lssc/L syn ~ Tp m[i/ c /z/b] ks 5tt for the typical synchrotron spectra with spectral 
index a = 1 between frequencies v\, in the infrared and v c in the ultraviolet. This does not say 
that the SSC mechanism is generally unimportant, it could still be the dominant process in the 
BL Lac objects where L 7 ks -L syn . 

2. 3 Spectrum 

All models do rather well in fitting spectral data (for an example, see Fig. 3). However, the 
sharp MeV break seen in some of the spectra is reproduced by EC and SSC models only by 
the ad hoc assumption of a break in the electron spectrum and it is difficult to find a physical 
process which could lead to this particular break. Coulomb losses must be negligible, otherwise 
the density of electrons is too large to comply with the weakness of the X-ray emission (absence 
of Sikora-bumro). Inhomogeneous EC mo delffl can easily reproduce MeV breaks assuming 
steep gradients along the jet. In the PIC model, the high-energy emission is due to unsaturated 
synchrotron cascades which, contrary to saturated cascades, can produce very hard X-ray spectra 
depending on the slope of the target photon spectrum. However, 7-ray attenuation is required 
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Figure 2: Redshift distribution of the 5 GHz to 100 MeV spectral index a rg for the strong EGRET blazars 
The 7/radio luminosity ratio is given by Lj/Lr, = (5 10 12 ) 1_QrE . 



to keep the 7-ray flux below Whipple limits. A natural source of the attenuating target photons 
is a heated dust torus. The PIC model predicts an average spectrum oc E~ 2 from MeV to TeV, 
and roug hly oc E~ 3 abovdEi The upper break is less certain, it depends on the range of distances 
in which proton cooling is important. In the PIC model, 7-ray emission is considered only from 
the part of the jet where the infrared target photons become optically thin. This is expected to 
be the most 7-ray luminous part of the jet. 



2. 



Cut-off 



Since the 7-ray spectra are broad featureless continua, the cut-off energy is of great diagnostic 
value. No cut-off is generally seen in the 7-ray spectra of blazars up to ~ 10 GeV above which the 
EGRET sensitivity drops rapidly. Generally, the flux level of the EGRET spectra extrapolated 
to air-Cerenkov energies exceeds the sensitivity of existing telescopes. Nevertheless, only 2(3) 
of the nearest blazars have been detected at TeV energies^ which is commonly ascribed to the 
attenuation of 7-rays duetp pair production in the extragalactic infrared background. The large 
zenith angle observations^ of Mrk421 indicate that the spectrum extends at least to 5-8 TeV 
with no evidence for attenuation. A possible excess of 7-ray showers at still higher energies from 
the stacked source positions of the nearest blazars (including Mrk421) with redshifts z < 0.07 
has been claimed on the basis of an advanced analysis of HEGRA array data. The significance 
of this excess has not yet been established unequivocally Eqs.(4) and (5) show a quite general 
result regarding the maximum possible energy in leptonic and hadronic emission models. If the 
observed 5-8 TeV emission from Mrk421 is to be explained by any model in which the primaries 
are accelerated electrons, the electron Lorentz factors must reach at least 10 7 . The same electrons 
are assumed to produce the radio-to-soft-X-ray spectrum by synchrotron emission. With an 
observed turnover energy of e ~ 300 eV and the formula for the characteristic synchrotron 
energy e ~ 10 -8 -B7 2 eV one obtains B ~ 3 10~ 4 G. Eq.(5) implies that the magnetic field in 
the acceleration region must be < 1 G adopting (optimistically) K = 0.1 and f c = 2 which is 
consistent with the above estimate. However, fields with a strength of B ~ 10 -4 G are found 
observationally in the hot spots of jets at the kiloparsec scaled. This is at least ~ 10 6 times 
larger than the size of the 7-ray zone inferred from variability. Adiabatic flux compression to the 
scale of the 7-ray emitting zone yields -B-fields somewhere in the range 0.1 — 100 G. Taking into 
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account a Doppler factor 5 ~ 10 does not change these conclusions significantly. Therefore, the 
magnetic field values implied by SSC models fitted to ~ 10 TeV emission seem unrealistically 
low. An estimate similar to Eq. (5) yielding ~ 10 TeV as the maximum possible energy for 
leptonic models has been calculated^ assuming ordered electric fields. The leptonic limit is in 
marked contrast to the hadronic models which naturally produce > 10 TeV energies. 



2. 5 Variability 

Doubling time scales of days to months seem to be common among strong EGRET sourceJl. 
Considering the long integration times required to obtain 5-a detections the minimum vari- 
ability time scales could also be shorter. Much shorter time scale variability has been discov- 
ered in Mrk421 at TeV energies taking advantage of the superior aperture of the air Cerenkov 
methoci. The shortest time scales expected for emission from blazar jets are of the order 
of At ~ rj/(c5) where the minimum jet radius can be estimated from asymptotic solutions 
of the equations for the radial structure of self-collimating hydromagnetic jet J]. This yields 
rj = lOIYs = 3 10 15 rio?7i8 cm where Tio = T/10 and nig = M/1O 8 M denote bulk Lorentz 
factor and the mass of the black hole, respectively, so that one obtains 

At min ~ ~ 5 10 3 m 8 s (9) 

CO 

assuming that blazars have favorably large Doppler factors 5 ~ T. In principle, one could further 
reduce this minimum time scale by another <5 _1 factor if the jet opening angle exceeds ~ 1/T or 
if the emission region is very thin and running through inhomogeneities. In spite of the short 
variability time scales, the pair creation optical depth of the emission region can be very small. 
The maximum optically thin energy for a spherical blob is given by E ~ 36f L^[e/6 2 ]At 4 TeV 
where L44 denotes the apparent luminosity in units of 10 erg/s and At 4 the variability time 
scale in units of 10 4 s. The low compactness is an essential requirement for the PIC model which 
would otherwise overproduce X-rays. In fact, the PICmodel has predicted emission above TeV 
from the presence of X-ray gaps in the blazar spectral. 
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Figure 4: Tentative score chart (model explains observations naturally [**], with ad hoc assumptions [*], not at 
all [0]). Blazar emission above ~ 10 TeV would practically rule out the leptonic models, as indicated by [(0)] 

2. 6 Correlations 

The claimed 7-ray/optical correlationil is consistent with EC models, if the variations are due 
to a varying bulk Lorentz factor. In this case, one would expect the amplitude of the 7-ray 
variations to be larger than the synchrotron variations, as seems to be the case in 3C27f3. SSC 
models predict proportional variations in this case and quadratic ones, if the variations are due 
to a variable comoving frame maximum Lorentz factor. The latter is suggested by UV/X-ray 
observations of blazars. The PIC model predicts basically the same variability pattern as the 
SSC model (protons instead of electrons scattering off synchrotron photons). The only difference 
is that variations of the electron maximum energy do not necessarily imply the same variations 
of the proton maximum energies, if the maximum energies are determined by a subtle balance 
between acceleration, energy losses, and expansion. Thus, the PIC model is less predictive in this 
respect. The inhomogeneous EC models, in which the time scales for 7-ray variations increase 
with energy and in which the TeV and optical emission sites are not in the same part of the jet, 
are clearly inconsistent with the claimed correlations. 

2.7 Warm dust ^y-ray attenuation 

As outlined in Sect. 1, the 7-ray emission models assuming external photons as the dominant 
target rely on a parametrization of the inner parsec of AGN as given by unified schemes. An 
element of the unified schemes which has more recently been recognized to be important is 
the dust torus heated by nuclear radiation feeding the accretion flow onto the black hole. The 
sublimation temperature of dust is ~ 1700 K and determines the inner edge of the torus to be 
at 

1 

r d ~210 4 p*Vr S (10) 

assuming a covering fraction of 0.5. This is roughly the same as the radius of the broad line 
region. The thermal NIR photons emitted from the warm dust torus absorb practically any 
7-ray with energy above ~ 300 GeV emitted from r < r<j. On the basis of a straw person's 
unified scheme, EC models are therefore ruled out as an explanation of TeV sources. 



3 Discussion and conclusions 



The results of Sect. 2 are summarized in Fig. 4. More refined models with more parameters 
could certainly improve some of the low scores, but some fundamental problems remain. The 
EC models have difficulties explaining the TeV emission due to the attenuating warm dust 
radiation field and the SSC models have difficulties explaining the enormous 7-ray luminosities 
of some quasars. The very attractive inhomogeneous EC models fail to explain 7-ray/optical 
correlations. The problems of the leptonic models could be ameliorated by assuming a transition 
from SSC to EC behavior when going from BL Lacertae objects to flat-spectrum radio quasars. 
Even if the emission properties can be explained by electron acceleration only, the absence of 
the Sikora X-ray bump is more in line with the energy transport in jets being dominated by 
protons rather than pairs and Poynting flux. 

In any model, the observed short variability time scales observed in Mrk42f push the re- 
quired acceleration times toward their limiting values allowed by diffusive shock acceleration 
and indicate that the size of the 7-ray emitting region is much smaller than its distance to the 
active nucleus. This is expected if the emission originates in the radiative downstream regions 
of shocks with a thickness of the order of the optical synchrotron cooling length. 

Observations above 10 TeV are crucial for the discrimination between leptonic and hadronic 
models, but they are unfortunately hampered by the expected cosmic attenuation of 7-rays at 
these energies even for nearby blazars. The detection of Mrk421 at 5 TeV and the possible 
cumulative excess of HEGRA events from the nearest blazars are encouraging hints that the 
cosmic transparency is not exceedingly strong. Ultimately, observations with high-energy neu- 
trino telescopes^ could bring us closer to an understanding of the puzzling nature of blazars 
and their relation to ultrahigh-energy cosmic rays. 
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